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Abstract 



>, 

^O , Anisotropies in the cosmic microwave background radiation due to gravity waves are in- 

Cn ■ vestigated. An initial spectrum of gravity waves may have been induced during an epoch of 

inflation. We study the propagation of such a spectrum in different Friedmann backgrounds 
with hot and cold dark matter, radiation and, possibly, a cosmological constant. We calculate 
its imprint as anisotropies on the cosmic microwave background, 
^i^ [ We also take into account that massless particles can source gravity waves by their anisotropic 

■~>^ ■ stresses. We consider general mixed dark matter models with and without cosmological con- 

1-^ ' stant. For a given, scale invariant input spectrum of gravity waves, we determine the de- 

pendence of the resulting spectrum of CMB anisotropies on the different parameters of the 



O ' model. 

H ■ 

'^ ' Keywords: Cosmology: cosmic microwave background; Gravitational waves; Cosmology: large- 

Ci ■ scale structure of the Universe. 



1 Introduction 

The theoretical and observational determination of anisotropies in the cosmic microwave back- 
ground radiation (CMB) has recently attracted a lot of attention. One has justified hopes to 
measure the CMB anisotropies to a precision of a few percent or better within the next ten years. 
Furthermore, if initial fluctuations are induced during a primordial inflationary period and no ex- 
ternal sources induce perturbations at later times, CMB anisotropies can be calculated by linear 
cosmological perturbation theory to very good accuracy. Since the detailed results depend not 
only on the primordial spectrum but also on the parameters of the cosmological model considered, 
the anisotropy spectrum may provide a mean to determine these parameters to an accuracy of a 
few percent. 

This problem has been extensively investigated mainly for scalar perturbations [|l|. 

As it is well known, also tensor perturbations can by generated during inflation. They play 
an important role: As shown in H, the presence of gravitational waves can crucially change the 



theoretical predictions of cluster abundance, which is an important test of cosniological models. 
In particular, power spectra of mixed dark matter (MDM) models normalized to the COBE 4- 
year data |Q provide cluster abundances higher than observed. This is one of the difficulties of 
standard MDM models. Taking into account a gravitational wave contribution, this inconsistency 
can be circumvented. Hence, the question how gravitational wave contributions depend on model 
parameters is very important. 

Here we discuss the model dependence of anisotropics due to gravitational waves for models with 
a total density parameter fi = 1 which are thus spatially flat. However, we vary the contributions 
of cold dark matter (CDM), hot dark matter (HDM) and a cosmological constant, which are given 
in terms of the parameters Qc i ^h and fl\. We also vary the number of degrees of freedom for 
massless neutrinos and hot particles. 

We consider a fixed input spectral index n = from inflation. For a general input spectrum 
{\h{tin, fc)p) = A(fc)^/c~^, our output spectrum (\h{t, fc)p) has to be multiphed by |^(fc)p. 

In Section 2 we present the perturbation equations and describe the models considered in 
this work. In Section 3 we discuss our results and we conclude in Section 4. The non-trivial 
relation between the temperature anisotropy spectrum, Ce, and the metric fluctuation spectrum 
(/ly (t, k)/i;m(t', k)*) for tensor perturbations is derived in the appendix. 

Notation: The Friedmann metric is given by a^{—dt^ + jijdx'^dx^), where a denotes the scale 
factor, t is conformal time, and 7 is the metric of a three space with constant curvature K. We shall 
consider a spatially fiat universe, the case K = exclusively. An over-dot stands for derivative 
with respect to conformal time t, while prime denotes the derivative with respect to kt = x. 



2 The models 



The basics of linear perturbations of Friedmann universes are discussed in H . We shall adopt the 
notation of pj in this work. We want to determine the evolution of tensor perturbations in spatially 
flat Universes which contain a fraction of cold dark matter (CDM), hot dark matter (HDM) and 
a cosmological constant A, such that: VLq = fin + ^c + ^A = 1, where fi, denotes the density 
parameter today, i.e. at Iq. We neglect the contribution of photons, massless neutrini and baryons 
(which may be included in CDM) to flo. 

The expansion of the Universe is described by the Friedmann equation for the scale factor: 

where pm is the total density of matter in the Universe, pm — Pc ^ Ph + Pva '^P-y- Here p^ denotes 
the density of photons, pi,g is the density of massless neutrini and pc, Ph denote the densities of 
CDM and HDM respectively 

The metric element for a Friedmann universe with tensor perturbations is given by: 

ds^ = a^{t){-dt^ + (% + 2hjj)dx'dx^) (2) 

where we choose c = 1, t is conformal time and a{t) denotes the scale factor. For tensor perturba- 
tions, the metric fluctuations hfj satisfy the conditions 

hf' = 0, hj'k' = (3) 

where k is the wave vector which may be set equal to (0, 0, k), such that the conditions (||) reduce 
to 

h{,^-h^^, and hl^hl^O. (4) 



We describe dynamics of tensor perturbations in a medium containing collision- less particles, 
whose anisotropic stresses are not damped by collisions. As long as the collision-less component 
is relativistic, it provides a source for gravitational waves. The evolution equation for tensor 
perturbations of the metric is given by g : 

hfj +2-hJj+ k^hfj ^ STrGa^pUij . (5) 

Here p is the pressure of the collision-less component and 11 denotes the tensor contribution to the 
anisotropic stresses, which in our case are due to the presence of relativistic, collision-less particles. 
Denoting the tensor part of the perturbed distribution function of the collision-less component by 
i<", n is given by 

1 /* 4 

n,j = ^ / —in^nJ - l/35,j)Fdvdn (6) 

where rii is a spatial unit vector, denoting the photon directions, v is the redshift corrected velocity 
and q the redshift corrected energy of the collision- less particles (see [5] ) . In the case of massless 
particles (massless neutrini) q = v. Liouville's equation leads to the following perturbation equation 
forFl, 

qF + vn-'kjF = qvn^rijh- ■'— — , (7) 

dv 

f denotes the unperturbed distribution function. 

The set of Eqs. (|l]) to (0) fully describes the evolution of tensor perturbations in media con- 
taining perfect fluids and collision-less particles. In our models we have in principle three kinds 
of collision-less particles: Hot dark matter, massless neutrini and, after recombination, the pho- 
tons. Studying the initial conditions, we shall find that for the growing mode anisotropic stresses 
are extremely small on super horizon scales. However, when the scales relevant for tensor CMB 
anisotropies (A ^ tdec) enter the horizon, t >> tdec HDM particles are already non relativistic. We 
may thus neglect their contribution to anisotropic stresses. Hence, we just consider the pressure 
anisotropy from massless neutrini and, after recombination, from the photons themselves. 

For massless particles we can simplify Eqs. (||) and (^ by introducing the brightness perturba- 
tion M 

M = / Fv^dv (8) 

Pi^ofl^ Jo 

In terms of M Liouville's equation (M) becomes: 

M + in^kjM = -An^n^hjj (9) 

and the anisotropic stresses are given by: 

Hy ^^ j {nrUj - l/3%)Mdr! . (10) 

Equations (p|) ; (Im) ; (W) for perturbations and eqn. (|^) for scale factor together with the usual equa- 
tions determining pc, pu Pvq and p~f form the closed system of ordinary differential equation which 
we have solved. 

We assume standard inflation according to which the initial amplitude of gravitational waves 
is independent of scale i.e., ( /i(ti„,A:)p) oc k~^ . 

Each solution of Eqs. (m),(M) and (M) can be presented as a sum of growing and decaying modes 
and an infinite number of modes corresponding to perturbations of the collision-less medium [pj . We 
are only interested in the growing mode which is given by the initial condition hjAt -^ 0) =:const. 
and h,j{t^ 0) = 0. 

If hij — 0, Eq. (^ does not admit a tensor contribution to M. In this case, M oc exp(m • k) 
and all components of the induced anisotropic stress normal to k vanish. Therefore, the correct 
(tensorial) initial condition for Af is M{t ^ 0) = and also Hy (i — > 0) = 0. 



These initial values remain unchanged as long as kt <^ 1. Assuming, e.g. spherical polarization 
and a flat spectrum from inflation, at some early time, kt <^ 1 for all wavelengths considered, we 
thus choose the initial conditions 

{\hl,\') ^ {\h{^?) ^ {\h\') = A^k-\ nn = ni2 = o, Af = 0, (11) 

where A is the amplitude of gravitational waves. It is easy see that on superhorizon scales (fct ^ 1), 
h =const. and the evolution of gravitational waves and as a result AT/T are independent of the 
model parameters. For scales kt fa 1, the metric perturbations begin to oscillate and eventually 
(kt 3> 1) damp away (see Figs. 1 and 2). The non-zero h then induces anisotropic stresses via 
Eqs.(^ and (10). Very often, these anisotropic stresses have been neglected in the literature. Here 



we find that their effect is indeed very small. There is typically about 1% additional damping due 
to the loss of some gravitational wave energy into anisotropic stresses. 

The main model dependence is the modification of the damping term (d/a) in the different 
backgrounds considered. 

We want to determine the CMB anisotropy spectrum induced by gravitational waves. Using 
that the brightness perturbation AI for photons is actually 

AT 

M = 4 . 

T 

we obtain by integrating Eq. (|9|) for photon;^ 

AT f*° ■ ■ 

—-(to, k, n) = exp(ik • nio) / n^n^h[j{t, k) exp(— ik • nt)dt . (12) 



The power spectrum, Ci of the CMB anisotropics can defined by the expansion of the correlation 
function into spherical harmonics. 

/AT AT ,A 1 



C(cos0) = (-^(io,a:o,n)--^(io,a:o,n')) = — E(2/ + l)QP,(cos0) (13) 



A somewhat lengthy calculation relates the gravitational wave spectrum \h{t, A:)p via Eqs. ( |l2|) and 
(O) to the C^'s. Taking into account the conditions (pi), it is possible to split this integral into two 
part coming from hj^ and /ifj- If we assume initially (|/i^iP) — (|^i2l) = 1-^1^ (no polarization), 
the two terms are equal. 

Ce = - f dkk^\I{£,k)\^ £{£-!){£ +!){£ + 2) , (14) 

with 

Ii£,k)^f%tHit,k) fl;^'^-J (15) 

(see alsojj), where je denotes the spherical Bessel function of order £. A self contained derivation 
of Eq. (^_5|) is presented in the appendix. 

Before we come to a description of the model dependence of the results, let us discuss the 
expected behavior in general. For a rough discussion we may neglect the anisotropic stresses in 
Eq. (ph. We first consider large scales which enter the horizon only after decoupling, kt^ec ^ 1- 
These scales contribute to \I{£,k)\'^ by roughly A^k~^j'j{kto)/£^. Inserting such a contribution in 
Eq. (nj) and integrating over k yields 

£2q ^ A^ . (16) 

^Prior to and during recombination, photons Thomson scatter with the electrons. For photons, we thus, in 
principle, have a collision term on the right hand side of Eq. (h). But this is only important on relatively small 

angular scales, i ~ 500, which we are not investigating here. We thus neglect the collision term. 



The integration over k is only justified if the main contribution to j| comes from the regime, 
where ktdec ^ 1, in other words, if kto ^ £ for a value of k with ktdec <C 1, which is equivalent to 

i < tdec/to -- 50. 

Perturbations on small scales, ktdec 3> 1 are damped by a factor of about (tenter/^rfec)^ ~ 
^/{ktdec)^ until decoupling^. Here tenter ^ 1/^ denotes the scale at which the mode k enters the 
horizon. Since the main contribution to Cg comes from the scales k with kto ~ £, we obtain an 
approximate behavior of 

fCi ^ ^2 f^\ for £ » tdec/h ~ 50 . (17) 

This analysis explains the generic behavior of the curves shown in Fig. 3. 

Let us now come to a more detailed analysis. Having calculated the metric perturbations hfj 
numerically, we can determine the CMB fluctuation spectrum according to Eq. ( [l5|) by means of 
numerical integration (we have used 60 to 100 point Gauss-Lagrange and Gauss-Laguerre integra- 
tions) and investigate it's dependence on the model parameters. 

3 Results 

We have solved Eqs. (J^), (|^) and dlOl ) numerically for tin < t < to choosing the initial conditions 
of the growing mode and unpolarized, isotropic waves, (|/ifiP) = (|/ii2p)- For a given model of 
inflationary initial perturbations, our results would have to be properly weighted and added to the 
scalar Q's. 

We have chosen a flat initial spectrum, such that hfn = Ak^^^^ and /ij„ ~ 0. 

We have investigated 80 models varying the flve parameters {ho,Q\,QH/^c,Piy, Ph), where 
ho is the Hubble parameter Ho in units of lOOkm/s/Mpc, /3^ denotes the number of degrees of 
freedom in massless neutrini and (3h is the corresponding number for hot dark matter particles. 

All the models lead to similar gravity wave induced anisotropics which, for reasonable parameter 
choices, differ by less than about 10%. The changes due to anisotropic stresses are extremely small, 
on the order of 1% or less. The main difference is caused by a non-zero cosmological constant, which 
enhances the damping at late times and thus leads to somewhat smaller perturbation amplitudes 
(see also [||). A similar effect is obtained if we increase the Hubble parameter. Hot dark matter 
does not induce significant changes since, at times when the wavelengths leading to substantial 
CMB anisotropies enter the horizon, hot dark matter is already non-relativistic, resulting in nearly 
the same expansion law as cold dark matter. In Fig. 1, we show the dependence of h{t,k) for 
fixed fc as a function of time varying several model parameters. We have chosen k — 20/to, which 
contributes mainly to an angle of ~ 6° in the sky, or to the C^'s with ^ ^^ 10 to 20. For some of 
these models we also show ft,(t, k) as a function of k for fixed time t = to/2 in Fig. 2. 

The solid line always shows standard CDM, i.e. Vtc = 1 and /3^ = 6 for comparison. The 
maximum amplitudes for standard CDM and CDM-f-A differ by more than 30%, while the mixed 
dark matter models show differences of about 1% only. The somewhat weaker damping due to the 
absence of the anisotropic stresses provided by a massless neutrino component and the decrease of 
a/a leads to the slight increase in amplitude for mixed dark matter models with ^h = 0.5. This 
result does not change if we increase the amount of hot dark matter. However, if we decrease VIh 
to 0.3 or less no amplitude change is left and only a small decrease in wavelengths builds up after 
about one oscillation. Due to the smallness of these effects, which remain of the order of 1% to 2% 
when translated into the C^'s, we disregard mixed dark matter models in what follows and claim 
that, on the level of 1% accuracy, MDM and CDM lead to the same gravitational wave spectra. 

^Here we neglect the short matter dominated period before decouphng and approximate the damping factor by 
its behavior in the radiation dominated era, {a/a) ~ 1/i. 



Changing the number of degrees of freedom in massless neutrini or HDM also induces very small 
differences of the order of 1%. 

Taking into account that an experiment always just measures the sum of tensor and scalar 
contributions and first has to disentangle the probably significantly smaller gravitational wave 
contribution, we can disregard such small effects, even if the experimental error is as small as 
possible, i.e. dominated by cosmic variance. 

The relevant parameters to be considered are thus ^Im = ^c + ^h, ^a and Hq. 

In the k dependence of h an additional effect comes into play: Due to the model dependence 
of io: the oscillations in ft, at a fixed fraction t = /tg of ig have different wavelengths if measured 
in units of ip ■ Models with a larger cosmological constant oscillate slower in Mq than models with 
small cosmological constant. Therefore, the cancelation due to oscillations in the integral (O) is 
more severe for models with small cosmological constant. This effect finally dominates over the 
larger amplitude of h which models with large A actually have. In Fig. 3 we show the C'i spectra 
for several models and in a Fig. 4 the relative differences are indicated. The A-models shown in 
frames (a) of Figs. 3 and 4 show a slightly increasing amplitude with increasing A. As can be seen 
in frames (b) of Figs. 3 and 4, increasing ho, which does not lead to an increase in the relative 
oscillation frequency, just decreases the fluctuations due to stronger damping. The detailed model 
parameters of the frames (b) are just given for information. The only parameters which really 
matter are VIa and ho. The variations induced by changing the other parameters are on the 1% 
level and thus swamped by cosmic variance. 

The variation of the tensor Cg spectrum for different cosmological models with fixed Hubble 
parameter which arc not already excluded by other observations than CMB anisotropics never 
exceed 10% for £ < 60, while variations of the Hubble parameter can lead to changes in the 
spectrum of up to 15%. 

4 Conclusions 

We have calculated the tensor contribution to the CMB anisotropics in mixed dark matter models 
with and without cosmological constant. We have included a previously neglected source term 
in the evolution equation for metric perturbations. Our findings are however quite modest: By 
reasons of cosmic variance, the statistical relative error in Ci measured from only one point in the 
universe is always Xj \J21 + 1. This is a very significant uncertainty, especially for the gravitational 
wave contribution which peaks around £ ~ 20 and has already dropped by a factor of about 2 at 
^ = 60 (see Fig. 3). 

In non of the considered models the influence of the anisotropic stress source becomes large 
enough to induce a difference in the Ci spectrum which is larger than cosmic variance. The same 
is true for hot dark matter contributions. Only an extremely large cosmological constant or a 
difference in the Hubble parameter can induce changes in the gravitational wave spectrum which 
are in principle observable but nevertheless small. 

This finding has one negative and one positive aspect: Unfortunately, the gravitational wave 
contribution does not contain detailed information about the cosmological parameters considered 
here and can thus not be used to measure them with high accuracy. On the other hand, since this 
contribution is so model independent, it conserves its information about the initial condition and 
thus about the amplitude and spectral index which it inherited during, e.g., an inflationary epoch. 
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Appendix: The C^'s from gravitational waves 

We consider metric perturbations which are produced by some isotropic random process (for ex- 
ample during inflation). After production, they evolve according to a deterministic equation of 
motion. The correlation functions of hij(\i.,t) have to be of the form 

(/ly(k,f)/l;;„(k,f')) = [hkjkikraHi{k,t,t') + 

{kikiSjm + kikmSji + kjkiSim + kjkmSii)H2{k, t, t') + 
kikj5imH3{k,t,t') + kikmSijH^{k,t',t) + 

+S^jSlrnHi{k,t,t') + iS^lSJm+S^.mSJl)H5{k,t,t')] . (Al) 

Here the functions Ha are functions of the modulus k = |k| only. Furthermore, all of them except 
H^ are hermitian in t and t' . This is the most general ansatz for an isotropic correlation tensor 
satisfying the symmetries required. To project out the tensorial part of this correlation tensor we 
act on hij it with the tensor projection operator, 

(A2) 

(A3) 



T, 



ab 



{PuP.m - (l/2)P„Pi™)P"«^P"' with 



P 

-t 7.1 






This yields 



Olmf^i'^j OilKjKjYi Oim'^l'^j ^jlf^ii^m Ojm'^l^i) ~r 






K rvj rv'i Kl K' 



(A4) 



From Eq. (12), we then obtain 



3 r /.to 

k^dkda 



/AT AT 

d^x — ^(n,x)— ^(n',x) 



V T 



T 



1 

2^ 



dt 



dt' exp(ik • n(to — t)) exp(— ik • n(to — t')) 



AT)* 



{hli'{t,-k)hyj*{t'Mn,n,n\n',^ 



(A5) 



/(k) 



d^A:exp(~ix- k)/(k) , 



Here dfij^ denotes the integral over directions in k space. We use the normalization of the Fourier 
transform 

If 1 

-= / d^x exp(ix • k)/(x) , /(x) = — — ^ 

where V is an (arbitrary) normalization volume. 

We now introduce the form ( A4) of < /i^-^^/i^-^' >. We further make use of the assumption that 
the perturbations have been created at some early epoch, e.g. during an inflationary phase, after 
which they evolved deterministically. The function Hc,(k^ t, t') is thus a product of the form 



H5{k,t,t') =H{k,t) ■H*{k,t') 



(A6) 



Introducing this in Eq. (A5) yields 



/AT AT 
( — (n) — (n) 



-!- j f k^dkd^^ [(n • n')^ - 1 + /i'^ +^i'^ - ^^i^Ji'{n ■ n') + ^V^] ' 
H[k, t)H* (fc, t') ex.p{ik^i{to - t)) exp(-ifc^'(io - t')) 



dt 



dt' 



(A7) 



where ^ = (n • k) and fi' — [n' • k). To proceed, we use the identity [^ 

exp((ifc^(io - t)) = ^(2r + \y jr{k{to - t))Pr{^i) 



(A8) 



r=0 



Here jr denotes the spherical Bessel function of order r and Pr is the Legendre polynomial of 
degree r. 



Furthermore, we replace each factor of /i in Eq. (A7) by a derivative of the exponential 
exp(i/c/i(to — t)) with respect to k{tQ — t) and correspondingly with //'. We then obtain 

r,r' 

2(n • n')^ I dtdt'jr{k{to ~ t))jr'{k{ta - t'))H{k,t)H*{k,t') 

dtdt'[jr{k(to - t))jr'{k{ta - t'))+j';{k{to - t))jr'{k{to - t')) + 

jr{k{to - t))f;,{k{t„ ~ t'))-j';{k{t„ - t))j';,{k{to - t'))]H{k,t)H*{k,t') 



4(n • n') / dtdt'jl{k{to - t))jl,{k{to - t'))H{k,t)H*{k,t') 



(A9) 



Here only the Legendre polynomials, Prilj) and Pr'{n') depend on the direction k. To perform the 
integration dVl^, we use the addition theorem for the spherical harmonics Yrs-, 



PM 



An 



- J2 Yr.{n)Y;,ik) 



(2r + l)^^ 
The orthogonality of the spherical harmonics then yields 

(2r + l)(2r' + l) /" df}^P^(^)P^,(^') 

r 
IdTT^Srr' J^ Yrs{n)Y;,{n') = 

s— — r 



(AlO) 



(All) 



In Eq. ( [A9| ) the integration over dfij^ then leads to terms of the form (n • n')Pr{n ■ n') and (n 
n')^Pr(n • n'). To reduce them, we use 



xPr{x) 



r + 1 
2r+l 



Pr+l 



2r + l 



-P 



r-l 



(A12) 



Applying this and its iteration for x'^Pr{x), we obtain 



AT AT* 
(^(n)^ (n')) 



-i^^(2r + l) f k^dk f dtdt'H{k,t)H*{k,t')( 



2(r + l)(r + 2) 
(2r- + l)(2r + 3) 



Pr+2 + 



-Pr 



2r(r - 1) 



(2r-l)(2r + 3) (2r-l)(2r + l) 



Pr- 



jriHto - t))jrik{ta - t')) - PrbriHta - t)j';{k{t(, - t')) 

+jr{k{to - t'))f;{k{to - 1)) - j';{k{to - t))f;,{k{to - m 

r + 1 



-Pr, 



-Pr-l 



]',{k{to~t))f,{k{to-t'))] 



(A13) 



_2r + l ^ 2r + l 

where the argument of the Legendre polynomials, n • n', has been suppressed. Using the relations 

r + 1 



3r = 



2r+l 



Jr+l + 



2r+l 



Jr-l 



(A14) 



for Bessel functions, and its iteration for j" , we can rewrite Eq. (A13) in terms of the Bessel 

functions jr-2 to jr+2- 

To proceed we use the definition of Ci : 



1 



/ AT AT \ 

(^(n)-^(n')) =— E(2; + l)C,P,(cos^ 



If we expand 



— -(n) = y^a^,my<;,m(n) 



(A15) 



(A16) 



and use the orthogonality of the spherical harmonics as well as the addition theorem, Eq. ( A10| ) 
we get 

We thus have to determine the correlators 



/AT*AT\ 



\ T T 



\r;„(n)y,,„,(n') 



(A17) 
(A18) 



Inserting our result (A13), we obtain the somewhat lengthy expression 



{aima.1 
2 



—SffiSmm' I dkk 

TT 



dtdt'H{k,t)H*{k,t')( 
Ji{k{to ~ t))ji{k{to - t')) X 



1 



(2£-l)(2£ + 3) (2^-l)(2^ + 3) (2£-l)2(2£ + 3)2 
4^3 4(^+1)3 



(2£-l)2(2£+l) (2<'+l)(2^ + 3)2 

[jiik{to - t))ji+2{k{to - t')) + ji+2{k{to - t))ji{k{to - t'))] X 

1 /2(^ + 2)(^+l)(2^2^2£-l) 2(^+ !)(<'+ 2) 8(£+ 1)2(^ + 2) 



2/ + 1V (2^-l)(2^ + 3)2 (2^ + 3) (2^ + 3)2 

[Mk{to - t))j£-2(fc(io - t')) + .n-2(k{to - t))je{k{to - t'))] x 
1 /2^(£-l)(2^2_^2£-l) 2^(^-1) 8i^{£-iy 



21 + 1 \ (21 -If (21 + 3) i2£-l)i2 [21 -\f 

+Je+2ik{to - t))Ji+2{k{to -t'))x 

f2{£ + 2){£+l) 4(£+l)(^ + 2)2 (£+1)2(^ + 2)2 



\{2£ + 1){2£ + 3) (2^+l)(2^ + 3)2 (2f + 1)2(2£ + 3)2 
+Je-2(k{to ~ t))je-2{k{to ~ t')) X 

2^(£-l) 4£{£-iy e^{£~lf 



(2£-l)(2£+l) (2^-1)2(2^+1) (2^-l)2(2£+l)2 



(A19) 



An analysis of the coefficient of each term reveals that this result is equivalent to Eq. (|lj) with. 

w„ , . ^ n+2{kito~t)) 2Mkito~t)) ji-2ikito~t)) 

^'' {2£ + 3){2£ + 1) {21 + 3){2e - 1) {2£ + l){2i - if ^ ^ 



jeikjto - 1)) 



(A21) 
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Figure 1 : The variable h is shown at fixed wave number k = 20/f o as funetion of time for different 
models. In frame (a), we consider models without HDM. The solid line shows pure CDM (with 3 
sorts of massless neutrini). The dotted, dashed, long dashed and dash-dotted lines show models 
with increasing Oa- In frame (b), mixed dark matter models with O// = 0.3 (dotted) and Q,h = 0.5 
(dashed) are compared with standard CDM (solid line). In frame (c) standard CDM (solid line), 
with Q.CDM = f^A = 0.5 (dotted line) and ^cdm = 0.5, JIa — fin — 0.25 (dashed line) are shown. 
In frame (d) we compare flA = 0.3 models with Hubble parameters ho — 0.5 (dotted) and ha = 0.75 
(solid); and J^a = 0.7 models with ho — 0.5 (dash-dotted) and ho = 0.75 (long-dash-dotted). 
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Figure 2: The variable h is shown at fixed time, t — to/2 as function of the wave number k for 
different models. In frame (a), we consider models without HDM. The solid line shows pure CDM 
(with 3 sorts of massless neutrini). The dotted, dashed, long dashed and dash-dotted lines show 
models with increasing J^a- In frame (b) standard CDM (solid line), with ^cdm = f^A = 0.5 
(dotted line) and ^cdm — 0.5, Vtjy — VLh — 0.25 (dashed line) are shown. 
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Figure 3: In frame (a), the angular power spectra of CMB anisotropies induced by gravitational 
waves are shown for models with different values for the cosmological constant. The solid line 
represents the model JIa = (solid line) and the amplitude increases with increasing A. In frame 
(b), we show the effect of increasing the Hubble parameter. The models chosen arc mixed dark 
matter models with cosmological constant, fln/^cMD = 0.15, fl\ ~ 0.7 with h^ — 0.5 (solid line) 
and /iQ ~ 0.75 (dotted line); and ^hI^cdm = 0.35, JIa — 0.5 with /iq = 0.5 (dashed line) and 
/lo = 0.75 (long dashes) are shown. 
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Figure 4: The relative differences for the models given in Fig. 3 are shown. In frame (a) the same 
line types as in Fig 3a are chosen, and the difference from standard CMB is indicated. In frame (b) 
The difference between /iq — 0.5 and Kq = 0.75 is shown for the model nn/^cDM ~ 0.35, J^a = 0.5 
(solid hne) and ^h/^cmd — 0.15, Q\ = 0.7 (dashed line). 
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